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Coherent one photon (2w) and two photon (w) electronic excitations are studied for graphene 
sheets and for carbon nanotubes using a long wavelength theory for the low energy electronic states. 
For graphene sheets we find that coherent superposition of these excitations produces a polar asym- 
metry in the momentum space distribution of the excited carriers with an angular dependence which 
depends on the relative polarization and phases of the incident fields. For semiconducting nanotubes 
we find a similar effect which depends on the square of the semiconducting gap, and we calculate 
its frequency dependence. We find that the third order nonlinearity which controls the direction 
of the photocurrent is robust for semiconducting tubes and vanishes in the continuum theory for 
conducting tubes. We calculate corrections to these results arising from higher order crystal field 
effects on the band structure and briefly discuss some applications of the theory. 

PACS: 42.50.Hz,42.65.Sf,61.48+c,78.20.Jq 



I. INTRODUCTION 



The magnitude and direction of photocurrents in semiconductors are ordinarily controlled using applied bias volt- 
ages. Interestingly the direction of a photocurrent in a semiconductor can also be controlled without bias voltages 
through phase coherent control of the incident optical fields. In a typical experiment an initial and final state are 
simultaneously coupled using two coherent excitations: one photon excitation at frequency 2u> and two photon exci- 
tation at to. The coherent superposition of these two excitations can lead to a polar asymmetry in the momentum 
space distribution of the excited photocarriers and therefore to a net photocurrent. The effect has been discussed 
theoretically JljH and observed experimentally in photoyield from atoms || and for photocurrents in semiconductors 
j|. Recently, two of us have proposed that for carbon nanotubes this effect could provide directional control of a 
photocurrent along the tube axis [^) and even suggests a novel method for biasing the diffusion of ionic species which 
intercalate within the nanotubes. 

In this paper we study the excitations which lead to this effect both for graphene sheets and for carbon nanotubes. 
In both these systems the low energy electronic properties relevant to most solid state effects are determined by 
an interesting feature of the band structure. The isolated graphene sheet has only an incipient Fermi surface; it is 
actually a zero gap semiconductors where the conduction and valence bands meet at precise points in momentum 
space. The carbon nanotube is a cylindrical tubule formed by wrapping a graphene sheet and for metallic tubes 
the "zero gap" feature manifests itself in a peculiar doubling of the low energy electronic spectrum with "pairs" of 
forward and backward moving excitations at both kp and —kp In either case the low energy electronic spectra 
are described by a two component Dirac Hamiltonian Q . 

In this paper we develop the theory of phase coherent one- and two- photon excitation within this model. The 
application to the graphene sheet turns out to be a useful pedagogical model which is unusual for a semiconductor and 
nicely illustrates the origin of phase coherent control of photocurrents for a graphene derived system. For graphene it 
is inappropriate to analyze the third order nonlinearity by analogy with the third order response in atomic systems, 
as has been done previously for semiconductors [Q. Instead we find that the third order response probes the rather 
unique geometry of the extended low energy electronic eigenstates which occur within the graphene sheet. The 
application of the model to a carbon nanotube shows, interestingly, that the third order nonlinearity is suppressed for 
excitations between the lowest subbands of any conducting nanotube and vanishes completely for transitions between 
the lowest subbands of a conducting "armchair" tube, but it is nonzero and robust for the gapped subbands of a 
semiconducting tube. In fact the effects we calculate are significantly stronger for semiconducting nanotubes than for 
a conventional semiconductor. In principle this effect might be used to distinguish conducting and semiconducting 
tubes in a compositionally mixed sample. Other possible applications of the idea will be discussed later in the paper. 

In this paper we briefly review the effective mass theory for the graphene sheet in Section II. In Sections III and IV we 
derive the interaction terms in the long wavelength theory which couple the electrons to time varying electromagnetic 
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fields and present a calculation of the coherent third order nonlinear optical excitations using this model. Section 
V applies the results to study the third order response of an isolated infinite graphene sheet. In Section VI we use 
the theory to study third order effects for conducting and semiconducting carbon nanotubes. A discussion and some 
applications of the results are presented in Section VII. 



II. EFFECTIVE MASS THEORY 

In this section we briefly review the effective mass description of the low energy electronic states. The theory is 
developed for an ideal graphene sheet, a section of which is shown in Figure 1. The primitive cell of this structure 
contains two atoms, labelled A and B in the figure. The lattice is unchanged after a translation by any combination 
of the two primitive translations vectors 



fi = o(l,0) 
T 2 = a(-,—) 



(1) 



where the bond length d = a/\/3. We introduce a pair of primitive translation vectors for the reciprocal lattice Gi 
such that Gi ■ Tj = 2^(5^ , yielding 
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(2) 



which generate a triangular lattice in reciprocal space. 
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FIG. 1. Direct and reciprocal space structures of the graphene lattice. The primitive cell contains two sublattice sites labelled 
A and B in the left panel. The right panel shows the first star of reciprocal lattice vectors and the first Brillouin zone. The long 
wavelength theory expands the electronic Hamiltonian for momenta near the K and K' points at the Brillouin zone corners. 

The critical points K(K') are important to our discussion, and they occur at the corners of the Brillouin zone of 
this reciprocal lattice at the positions 



K 



\{Gi + 2G 2 ) 
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K' = I(-G 1+ G 2 ) 
4tt 1 V3 



(3) 



The "bonding" and "antibonding" tt electron bands meet precisely at these K(K') points in reciprocal space. This 
band touching is required by symmetry for this system and it is correctly described by the simplest model for electron 
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propagation within the graphene sheet which is a tight binding model in which the hopping of an electron between 
neighboring sites is set by a single energy, t. Thus we have 

hfa/ = (^l-ff = t (nearestneighbor /iu) 

(otherwise) (4) 

Working in the sublattice basis and at crystal momentum k we have the Hamiltonian 

/ n l + e -tk-T 2 + -ik-(f 2 -T 1 )\ 

If we set k = K + q and expand the Hamiltonian for ga << 1 wc obtain the long wavelength Hamiltonian 

y/3ta ( q x + iq y 

= TivpS* ■ q (6) 

where a are the 2x2 Pauli matrices. A similar expansion near the K' point yields 

Hk> — —fivpd ■ q (7) 

Identifying each of the critical points with the index a so that a = 1 denotes the K point and a = — 1 denotes the K' 
point, these Hamiltonians can be rotated into diagonal form with the unitary operators 

Ua ® = 71 ( -ae~^ ae~^ ) (8) 

where 9 = tan _1 (g y /q , x ). Thus 

Ul{q)H a {q)U a {q) = hv F ( ~^ ^ ) (9) 

so that equations (6) and (7) describe pairs of bands which disperse linearly away from the critical K and K' points. 
Note also that H a (q) = H^(—q) as expected. Equations (6) and (7) and the unitary rotations in equation (8) provide 
the appropriate description of all the low energy electronic excitations required for this problem. 



III. GRAPHENE-FIELD INTERACTION HAMILTONIAN 



In this section we collect several results we need to describe the coupling of electrons described by Equations (6) 
and (7) to the electromagnetic potentials. For a particle of charge Q interacting with the electromagnetic vector 
potential A and scalar potential $ the momentum and energy are shifted p — > p — QA/c and E — > E — Q$ 0j. Thus 
the interaction Hamiltonian which couples the Dirac particle to the time varying vector potential (A x , A y ) is 

H a< i n t = -av F Q(A x a x - aA y a y )/c (10) 

Note that this interaction operator can also be obtained by calculating the velocity operator as the commutator of 
the position operator with the unperturbed Hamiltonian: 

v a =^[H a ,f^\ (11) 

and therefore 

v a = av F {(7x, —ot<7y) (12) 

so that H a j nt = —ja ■ A/c. 

It is useful to rotate the interaction Hamiltonian (10) into the band basis using the unitary operators in equation 
(8). To do this we write q = (q x ,q y )/\q\ and compute H b a mt = W a {q)H a ^ nt (q)U a {q) giving 
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c \ taz ■ (qx A) q - A J 

This demonstrates that the coupling between the Bloch electrons and the vector potential depends on the angle 
between v and A and that the intcrband matrix elements (which are the off diagonal terms in Equation (13)) vanish 
when the two are collinear. Indeed, H a (q) and H a ^ nt {q, A) commute along these special lines in reciprocal space, so 
that interband transitions are forbidden along this trajectory. This peculiar feature can be traced to the absence of a 
mass term in the effective Hamiltonians in equations (6) and (7) which would ordinarily mix the plane wave solutions 
to (6) and (7) and thereby allow interband transitions by coupling with the long wavelength current operator. When 
v and A arc not collinear intcrband transitions are allowed and the transition amplitudes are fixed by the mismatch in 
their orientations in the graphene plane. This has interesting consequences for possible coherent control of nonlinear 
optical processes in the nanotubes, as we show below. 



IV. NONLINEAR OPTICAL EXCITATIONS 



In this section we present a calculation of the transition probabilities for the third order nonlinear optical excitations 
among the electronic states given by the models in Sections II and III. We introduce time varying fields of the form 

A(r, t) = A LU e- luJt+t ' pl + A 2LO e- 2wt+l ^ + c.c. (14) 

and study the response of the system to third order in these exciting fields. Asymmetries in the photocurrent are 
controlled by the coherent excitation of electrons from an initial state to a final state by one photon (2u>) and by two 
one photon (u>) processes. The coherent mixing of these two processes is studied by evolving the density matrix to 
third order in the exciting fields and isolating the terms proportional to A^A^A^u- 

It is convenient to study the time evolution of the one particle density matrix p — (\&t(r)\& (r')). The Hamiltonian 
for our system is H a + H a ^ nt and we work in the Heisenberg representation so that 

f t =\[H a , int {t),p] (15) 

where H a is the free particle Hamiltonian and H a . int (t) = e lHat H a ^ nt e~ lHat . In the band basis the density matrix 
in the initial state has the form 

Po = P(t = -oo) = (J °) (16) 

since only the negative energy states of the Hamiltonians (6) and (7) are initially occupied. 
We expand the density matrix order by order in the exciting fields 

p(t) = Po + Pl + p 2 + P3 + - (17) 

Integrating equation (17) to first order in the applied fields gives 

iev F ( Fi(t) 
— ' 

where 



* = liT I -*T(t) o ) ( 18 > 



az ■ {q x A u] )e- l ( A+ ^ t + l <t> 1 az • (? x l^e-'^)'- 
= a Tx + 



— iS —A + u> — iS 

az-(qx l 2 ^) e -'(A+2^)w<fe gz.(qx A^)e^ l(A ^ 2u))t ^ 
+ -A -2uj — id + -A + 2lu - id ( ' 

where A = 2vpq and S is a positive infinitesimal. 

The second order terms pi include the lowest order changes to the occupation probabilities which can be induced 
with excitation by the u or the 2uj fields 



2lTP 2 V 2 

(h)22 = -{h)n - —r^{\ocz ■ (q x A U )\ 2 S(A -w)+ \az ■ (q x A 2uJ )\ 2 S(A - 2u>)) (20) 
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as well as oscillating nonlinear off diagonal coherence terms 



2„,2 



( P 2)i2 = - l j^F 2 (t) (21) 
Anticipating the situation 2w « A, the most important contribution to -F^i) has the form 

2U {-A + u>-iS)(-A + 2u>-iS) [ ' 

The second order coherence term in equation (22) leads to a transition rate which is third order in the applied fields 
and is the source of the polar asymmetry of the photocurrcnt. 

(h)22 = Re ((^ ' « x ^u,)(« ' ■ Q x A^)e^- 2 ^)S(A - 2cu) (23) 

An c? 

Equation (23) contains a factor of two from the sum over the (physical) spins. Equation (23) presents the main result 
of the paper. It shows that the transition rate depends on the polarization and phases of both exciting fields and 
the Bloch wavevector q. We will explore the consequences of the geometric structure of this result for the graphene 
sheet and for carbon nanotubes in the following two sections. For the moment we note that the result is odd in 
the direction of the Bloch wavevector q and even in the critical point index a (it depends on a 2 ) and therefore the 
symmetry breaking nonlinearity is nonzero after integration over the full Brillouin zone. 



V. APPLICATION TO GRAPHENE 



In this section we apply the formalism developed in section III to study the coherent optical control of photocurrcnts 
for a single graphene sheet. The model nicely illustrates the selection rules which apply in this geometry, and the 
results can be extended to analyze the more complex situation for the nanotube which will be presented in section 
VI. 

We note that both equations (20) and (23) contain terms which describe transitions from the valence to the con- 
duction band at the frequency 2lo — A. Equation (20) is the ordinary linear absorption in the material. Interestingly 
we see that the angular distribution of the excited photocarriers is not isotropic but rather follows a sin 2 </> dependence 
with respect to the polarization of the exciting radiation. Nevertheless this angular distribution has even parity and 
thus does not produce a net current. On the other hand equation (23) gives an angular distribution that breaks the 
inversion symmetry of the graphene sheet. The symmetry breaking is actually implicit in the coherent superposition 
of the exciting fields. We will estimate the prefactors to compare the relative strengths of these terms for accessible 
laboratory fields later in the paper; for the moment we note that the nonlinear terms in equation (23) typically con- 
tribute ss 10~ 3 of the total transition rate, and thus the induced anisotropy while nonzero (and we believe measurable) 
is a subtle effect. 

The angular distribution in equation (23) is controlled by the relative polarizations and phases of the incident 
and A2ui fields. Figure 2 displays polar plots of the angular distributions for the situation where the amplitudes 
are in phase (i.e. A^ and A^ u in equation (14) are presumed to be real) for various incident polarizations. Note 
that the underlying electronic dispersion relations are completely isotropic in the linearized theory, and thus only 
the relative polarization of the two exciting fields is relevant for the interference pattern. In all cases A^ is taken to 
be polarized along the horizontal direction shown by the arrow in the plots. In each plot we observe a node in the 
current distribution along this direction. This follows from the symmetry of int in equation (13) which shows that 

interband coupling is prohibited for q parallel to A. Nonetheless, the situation for collinear u and 2uu excitation clearly 
shows the asymmetry between the "forward" and "backward" distribution of the photocurrent. The situation is more 
interesting when the exciting fields arc noncollincar. We observe that the angular distribution develops a "three" lobe 
structure. Ultimately when the exciting fields area mutually orthogonal, we recover the "two lobe" pattern with the 
angular distribution rotated by 7r/2 with respect to the polarization of the incident u field. It is useful to quantify 
the anisotropy of the distribution by calculating the average polarization of the net photocurrent (cos <j>) and (sin <j>) 
averaged over this distribution. One finds 

(cos <f>) = ^ cos 6 

(sin </>} = ~ sin6 (24) 
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so that when the 2uj field is tipped by an angle 9 with respect to the u field, the photocurrent is oriented in the 
direction —9. Finally, the "sign" of the effect is determined by the relative phases of the two exciting fields. Note 
that the phase delays <j>\ and 02 in the exciting fields of equation (14) modulate the transitions rates j| in equation 
(23) in the combination 

(P3)22 -* (^3)22 cos(0 2 - 20i ) (25) 

This does not change the qualitative features of the angular distribution but it can modify both its magnitude and 
its sign. Thus the angular distribution in the collinear case A9 — can be reversed by advancing the phase of the u> 
fields by 7r/2. 




A8 = A9 =7C/6 A8 =7C/4 A9 =n/3 A8 =71/2 

FIG. 2. Angular distributions of the transition rates given by equation (23). In each panel the A u is polarized along the 
horizontal direction (the direction of the arrow in each plot) and A9 is the angle between the A^u and Au, fields. The polar 
plot gives the transition rate as a function of the angle of the Bloch wavevector q with respect to the direction of A^. These 
angular distributions are superimposed on the angular distribution for the direct transition rate which is given by the second 
term in equation (20). The solid curves correspond to (403)22 > 0, dashed curves to (^3)22 < when cos(02 — 20i) > 0. 



VI. APPLICATION TO NANOTUBES 
A. Low Energy Theory for Semiconducting Tubes 

In this section we apply the formalism of section III to study phase coherent control of a photocurrent on a 
carbon nanotube. The essential complication is that the wrapped structure of the nanotubc quantizes the allowed 
crystal momenta so that the transition rate automatically contain an intrinsic anisotropy. Nevertheless, the formalism 
developed in section III can be extended to include this situation. 

We first define the geometry for the single wall nanotube. The nanotube is a cylinder formed by wrapping graphene 
sheet and the wrapping can be defined by the graphene superlattice translation vector around the tube waist. We 
adopt the primitive vectors of equation (1) as a basis and represent the superlattice translation Tmn as 

T M N=MT 1 +NT 2 = (M + j,^^)a (26) 

It is useful to define two unit vectors defining the longitudinal and azimuthal directions within the graphene plane 

e; = (cos# c ,sin# c ) 

e a = ( — sin# c , cos C ) (27) 



where 9 C = cos" _1 (M + N/2) /(VM 2 + N 2 + MN) is the chiral angle of the tube. The wrapping of the tube quantizes 
the allowed momenta along the azimuthal direction k ■ e a — 2-nnj '(av 'M 2 + N 2 + MN) while the electrons obey free 
particle boundary conditions along the tube direction and the longitudinal component k ■ e; can take any value p|-^0[ . 
Thus the loci of allowed momenta are "lines" in reciprocal space. These lines need not intersect the critical K and 
K' points which are used as a reference for the long wavelength theory. To determine the mismatch between the 
allowed crystal momenta and the K and K' point wavefunctions we resolve the Bloch wavevector at K and K' into its 
longitudinal and azimuthal components. We find 

3a 2' 2 ' V ' 
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and 



3a \2VM 2 + N 2 + MNJ V ' 

which lies along the locus of allowed wavectors when 2aM + (3 + a)N = 6n. One third of the (M, N) tubes satisfy 
this condition, and for the remaining two-thirds of the tubes the K(K') momenta are mismatched to the kincmatically 
allowed momenta by a mimimum amount 

^_ j^mod(2aM+(3+a)iV,3) ^q-j 



3aVM 2 +N 2 + MN 

Representing the "reduced" Bloch wavector with the complex number q = q x + iq y = qe l8a and the momentum 
mismatch by A a = iAe l9c after a rotation of the coordinate system by the chiral angle 9 C (so that the x-axis runs 
parallel to the the tube length) the Hamiltonians in equation (6) and (7) can be written 



Note that in equation (31) H a (q) = H*L a (—q). The spectrum is now E(q) = ±y/q 2 + A 2 and the Hamiltonian is 
diagonalized with the unitary transformation 

U °M ^ 71 ( - J"- J«n ) ( 32 ) 

where 7 = tan _1 (A/g). This is the rotation identified in Equation (8) for the unfolded graphene sheet with A playing 
the role of the y-component of the momentum. With this identification the interaction Hamiltonian for the nanotube 
analogous to Equation (13) in the band basis is 

h ^ _ ev F A 1 ( —q iaA 



<- ( ^ ) ^7^tM--A q ) ^ 

Note that the off-diagonal terms which describe the amplitudes for interband transitions depend explicitly on the size 
of the semiconducting backscattering gap A and vanish for the lowest subbands of a conducting nanotube as shown 
in Figure 3. 




FIG. 3. Optical excitation between lowest subbands of a conducting tube (left panel) are forbidden in the long wavelength 
theory. They are allowed for a semiconducting tube (or for the gapped subbands of a conducting tube) as shown on the right. A 
denotes the crystal momentum mismatch between the valence and conduction band states and the K(K') points of the graphene 
sheet. The right hand panel illustrates one- and two- photon excitations which interfere to produce a polar asymmetry in the 
photocurrent. 

Thus when the exciting fields are collinear and directed along the tube direction the third order transition rate is 
(Ps)22 - " - / A 2 qRe(A 2oJ A^A^e^- 2 ^)5(2E(q) - 2hu>) (34) 
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Equation (34) is the origin of the asymmetry discussed in reference || . We note that the result depends on the square 
of the magnitude of the gap A and it vanishes for transitions between the lowest subbands of a conducting tube. 
The result is odd in the reduced momentum q which produces the asymmetry between forward and backward moving 
photocarriers. The third order transition rate is very small for high exciting frequency since the high energy electrons 
are very weakly backscattered through the mass term in equation (31) and behave essentially as free particles. 
These properties are displayed in Figure 4 which shows the third order transition rate between two bands of a 
semiconducting tube as a function of the exciting frequency. It is interesting to note that the expected divergence 
in the one dimensional density of states at thrcshhold is exactly cancelled by the momentum prefactor q in equation 
(36) and thus the spectrum shows only a step-like singularity at the threshhold. 




500 1000 1500 2000 
fico (meV) 



FIG. 4. Frequency dependence of the third order transition rate leading to anisotropy of the photocurrent. The solid curve 
gives the transition rate of Equation (30) as a function of the exciting frequency hui. To display the spectra we have taken a 
semiconducting gap ?ivfA — 800meV and normalized the incident intensity so that \A\ — 10~ 9 T-m at all frequencies. The 
dashed curve is the transition rate weighted by the final state velocity. Band edge states are strongly scattered by the mass 
term and do not contribute effectively to the photocurrent. 

Thus the transition probability for right and left moving photocarriers jumps discontinuously across the critical 
point at q = 0. Note however that the states near the bandgap have no group velocity and cannot contribute to 
the photocurrent so the velocity weighted transition rate (which is more relevant to this application) goes to zero 
at threshhold. This is shown by the dashed curve in Figure 4. The results in Figure 4 show the nonlinear injection 
rate for a perfect defect free semiconducting tube. A slowly varying impurty potential (long range disorder in the 
notation of reference |]l2|| ) can produce an additional channel for backscattering and will therefore further suppress the 
group velocity for electronic states near the band edges. In the presence of disorder we therefore expect an additional 
rounding of the current injection rate, similar to that shown as the dashed curve in Figure 4. The range and strength 
of this suppression will depend sensitively on the details of the long range impurity potential. The results of Figure 
4 are insensitive to this additional backscattering deeper into the particle-hole continuum. 



B. Crystal Field Effects For Conducting Tubes 

Equation (34) gives the third order nonlinear response in the long wavelength limit where we can linearize the 
electronic bands around the critical K(K') points. Corrections to this result can be obtained in an expansion in 
qa and physically arise from crystal field ("trigonal warping") effects in the underlying band structure. The most 
significant such corrections occur for conducting tubes. Equation (36) gives a vanishing transition rate for excitations 
between the lowest bands of a conducting tube and trigonal warping of the band structure of the graphene sheet 
provides a mechanism to "turn on" these transitions even for conducting tubes. Thus a third order nonlinear response 
is symmetry allowed for the lowest subbands of a conducting nanotube, though it strictly vanishes in the long 
wavelength limit we have discussed so far. 

To investigate the trigonal warping effects we re-derive the interaction Hamiltonian without adopting the effective 
mass representation. To do this we note that in the presence of a vector potential A the Hamiltonian (5) is perturbed 
through the Peierls substitution k — ► fc — Q-A. Therefore we can calculate the current operator using = —dH/dA^ — 
Q-dH/dk^. In the site representation the Hamiltonian has only off diagonal elements, so we can write 

H int = - e -i.( ° vm) (35) 
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We also note that the Hamiltonian is diagonalized with the unitary transformation: 



£/(£) = i=( t . (jE) j (36) 



V2 v 



|t(fc)| |t(fc)| 



which is the discrete lattice analog of the continuum result in equation (8). Thus we can rotate the interaction 
Hamiltonian into the band basis according to W (k)Hi nt (k)U(k) which gives 

H b e A ^ f -*te(t*d ku t) ilm(t*d k J)\ 

mtl j ~ c\t(k)\ \-Um(td k)l t*) Re(t*d K t) J [6 '> 

Explicit evaluation of the matrix elements in equation (36) for a general chiral nanotube is complicated. In general 
one may have interband matrix elements between lowest subbands of a conducting tube (which are the off diagonal 
terms in equation (37)); albeit with greatly reduced magnitudes - the scale of these matrix elements is typically 
~ 10 -2 the scale for the matrix elements in (31) which are produced by the mass term in the linearized theory for a 
semiconducting tube. An important exception to this rule for conducting tubes occurs for the armchair (M,M) tubes. 
Then one finds that t(k) = e 27Tl ^ 3 (l + 2cos(k x a)) for propagation in the lowest subbands of the tube, and we have 

{t*(k x )d k t(k x ))/\t(k x )\ = iasgn(l + 2cos(fc x a))sin(fe c a) (38) 

Thus near the critical points the diagonal elements of the velocity operator are ±vf and the off diagonal components 
vanish everywhere. Note that this occurs because of a tube symmetry; the armchair tube retains a mirror plane which 
contains the tube axis so that the two lowest subbands of the conducting tube can be indexed as even or odd under 
reflection through this mirror plane. The vector potential along the tube axis is even under the mirror reflection and 
cannot couple even and odd subbands. On the other hand, for a zigzag tube one has t(k) = e lk v a /^{\ -\- e -«v / 3fc H a/2-j 
This vanishes for k y = 27r/\/3a which corresponds to a "face" of the Brillouin zone in Figure 1. Thus one finds for 
the conducting zigzag tube: 

ta (4= cos 2 (\/3fca/4)) - # cos(V3ka/4)e^ ik y a / 4 ) 

(t*d k t)/\t\ = 22 2 (39) 

V >n ' |cos(V3fca/4)| V ' 

Therefore near the crossing point k y = (2n/y/3a) + q one finds 

(40) 

Thus the diagonal matrix elements of the velocity operator (the real part of equation (39)) are constant (vF(l+0(qa) 2 ) 
while the off diagonal elements (the imaginary part) vanish proportional to qa near the Fermi points. This implies 
that the product of the matrix elements in the third order transition rate vanish as (qa) 2 for the conducting zigzag 
tube. This changes both the magnitude and the frequency dependence of the third order transion rate. We obtain 

3 

(p 3 ) 22 = ^^a 2 UB^(q)Re(A 2u A- u A^e i ^-^)S(2(\t(q)\ - M) (41) 

The result is plotted in Figure 5 using the same normalization as in Figure 4 for comparison (note the scale change). 
One finds that the transition rate vanishes linearly in frequency, and is suppressed by f=a 10 -2 with respect to the 
interband transition rate for a semiconducting tube. This reflects the fact that at low energy the effects of trigonal 
warping are relatively small compared to the backscattering from the mass term in the low energy Hamiltonian for 
a semiconducting tube. We note that calculations of the frequency dependence of the resonant Raman cross section 
for conducting tubes show a strong enhancement of the cross section near the first interband threshhold, also 
demonstrating the suppression of interband transition matrix elements between the lowest conducting subbands in 
these structures. 
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' 500 1000 1500 2000 

■firo (meV) 

FIG. 5. Frequency dependence of the third order transition rate involving the lowest subbands of a conducting zigzag tube, 
which produce an to anisotropy of the photocurrent. The dashed curve uses a linear dispersion relation for the electronic states 
with the matrix elements are computed using the lattice theory of equation (37). The normalization of the incident fields is 
the same as for the results of Figure 4 so that the rates can be directly compared (note the vertical scale change ). These 
interband excitations are symmetry forbidden in the Dirac theory but become weakly allowed in the presence of crystal field 
effects on the low energy electronic states. For a conducting tubes third order transitions between the gapped subbands (not 
shown) provide a much stronger nonlinear third order response, as shown in Figure 4. 

Figure 5 presents only the results for excitations coupling the lowest subbands of a conducting tube. Transition 
rates between gapped subbands are described by Equation (34) so that the transition rate displayed in Figure 4 should 
be superposed on these results. This situation calculated for a zigzag tube illustrates the generic behavior for a a 
general (M,N) tube if one wishes to calculate beyond the linearized theory. Analogous results for arbitrary chiral 
tubes can be obtained by direct evaluation of the matrix elements in equation (37). 

C. Noncollinear Fields 

This treatment can be extended to include the situation where the exciting fields are not collinear. Interestingly 
this does not change the qualitative frequency dependence shown in Figure 4, although the prefactor is altered for 
noncollinear fields. 

We will consider only the case where the exciting fields are orthogonal, since any incident field can be resolved into 
its longitudinal (along the tube) and transverse (perpendicular to the tube) components. We observe that for a field 
perpendicular to the tube axis we have allowed interband transitions only when the azimuthal quantum number m 
changes by ±1 since the vector potential A "seen" in the tangent plane of the tube is A ■ <p — Acos<p where <$> is 
a unit vector which circulates counterclockwise around the tube waist. For the graphene sheet this is equivalent to 
introducing a spatially varying vector potential with wavevector 1/R where R is the tube radius. Thus the third order 
nonlinear process we are seeking is symmetry forbidden if A 2uJ has a transverse polarization (the lowest subbands 
have the same azimuthal quantum numbers.) However, it is possible to have the situation shown in Figure 6, where 
A w is perpendicular to the tube axis, and A 2( ^ is polarized along the tube direction. Here the virtual intermediate 
state for the two-photon process is provided by a higher azimuthal subband. 
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FIG. 6. Intefering excitations when the u) field is polarized perpendicular to the tube and the 2ui field is polarized along the 
tube direction. The 2u field excites transitions between subbands with the same azimuthal quantum numbers. The uj field 
excites transitions with 8m = ±1. 



This reduces the strength of the effect, but not the overall frequency dependence which is controlled by dispersion 
of the lowest azimuthal subband that is accessed to second order in A u . 

We modify the interaction Hamiltonian equation (33) for the situation where the exciting radiation is polarized 
perpendicular to the tube direction. In the "site" basis one finds that the interaction Hamiltonian for this polarization 
is 

ev F A ( — i 



where A = A(y) = AqCos(^). The y dependence implies that ths interaction couples subbands with a difference in 
azimuthal quantum numbers m such that 8m = ±1 and we will explicitly condsider only the two low energy pairs of 
subbands as shown in Figure 6, which we label 1 and 2. The Hamiltonian in equation (42) can now be rotated into 
the band basis using the unitary rotations of equation (32) in the combination H h a int — ul(q)H^ int U\{q) which gives 

b _ evpA i(72 _ 7l)/2 / sin(a( 7 i + 72 )/2) i cos(a( 7 i + 72 )/2) \ 
a ' mt c \ -icos(a(7i + 72 )/2) sin(a( 7 i + 72 )/2) J 



Thus the second order coherence term (analogous to equation (22) for the graphene sheet) is 
(p 2 ) 12 = > 2 ^ a 2 sin( 7l+72)e -^-^ 



The coherence factor appearing in (44) can be re-expressed in terms of the Hamiltonian parameters 

S in( 7l+72 ) = ^V F M ( 45 ) 



where E m (q) = ^ 'q 2 + m 2 A 2 and A is the lowest gap of the semiconducting tube. Note that the effect vanishes for 
Ai = A 2 i.e. between subbands of the same azimuthal symmetry. The second order coherence factor leads to the 
symmetry breaking third order transition rate 

(Pa) 22 = 4n fff (A2 ~ A F l)Al9£l Re(^A_^_,e^-^ 1) ^(2(E (g) - M) (46) 
which has exactly the same frequency dependence as the result of equation (34). 



VII. DISCUSSION 



Third order phase coherent control of photocurrents have been studied and demonstrated for semiconductors (e.g. 
GaAs ||Q) and since the effects calculated for carbon nanotubes are strongest for semiconducting tubes, it is 
appropriate to compare these effects. We find that the predicted effects are significantly stronger in nanotubes than 
for conventional semiconductors. This occurs because of the larger carrier velocities and the longer carrier relaxation 
times which are expected for the nanotubes. For carbon nanotubes this is particularly interesting since this third 
order nonlinearity provides a method for current injection without contacts. It has proven experimentally difficult to 
fabricate low resistance electrical contacts with carbon nanotubes by conventional submicron lithographic methods. 

For an incident intensity S — 10 2 kW/cm 2 the electric field amplitude E = 8.5 x 10 5 V/m « 10 6 V/m. At an 
optical frequency u> — 10 15 s -1 this corresponds to a vector potential amplitude \A\ « 10 9 T — m (which is the value 
used in obtaining the results in Figures 4 and 5.) Then we find that the typical carrier injection rate dn/dt sa 10 6 s~ x 
per unit cell. For hot photo-excited carriers the relaxation rate is presumed to be dominated by phonon emission, 
for which we estimate a carrier relaxation time r ss 1 ps so the steady state distribution gives n ~ 10~ 6 carriers 
per tube unit cell (note that the unit cell contains typically 40-60 carbon atoms around the tube circumference). 
Summing over the electron and hole contributions to the photcurrent and over the two electronic branches (K and 
K') wc obtain an induced current I » 0.4 nA, or an effective current density J ~ 260/^A//mi 2 . This is 10 — 10 2 



11 



larger than the induced density predicted for third order transitions between the valence and conduction bands in 
GaAs The enhancement is due mainly to the relatively large carrier velocity for the carbon nanotubes, and 

the larger estimated carrier relaxation times. For conducting tubes this enhancement is partially offset by the small 
interband matrix elements between the lowest subbands of conducting tubes; for a conducting zigzag tube we estimate 
the photocurrent density J s» 5/iA//im 2 under the same assumptions, a value which is comparable to that found for 
conventional semiconductors [jl3|. We note that even with these long carrier relaxation times, one should be able to 
achieve a steady state distribution during a 100 ns incident pulse. For conducting nanotubes one has the additional 
difficulty of resolving this signal over a background free carrier density fit, « 10~ 5 produced by ordinary one-photon 
excitation between the lowest subbands (first term in equation (20)). Since this is a "non-polar" contribution, i.e. it 
does not contribute to the photocurrent, the nonlinear contribution can be identified, in principle. 

The angular distributions calculated for interband excitations in graphene sheets show a similar structure to the 
angular dependence calculated for the third order rate for transitions from the heavy hole band to the conduction 
band in GaAs (h|] . In both cases the net induced current is polarized along the direction of the exciting field, but 
the current distribution is peaked away from the field direction. The high symmetry of the graphene sheet provides 
an additional interesting degree of freedom, namely control of the direction of the injected current by controlling 
the relative polarizations of the incident fields, as displayed in Figure 2. It would be very interesting to carry out 
experiments on graphite (either bulk or thin films) to verify the predicted angular dependence. Quantitative studies 
of the magnitude of the effect would be very useful as a probe of the scattering processes which control the dynamics 
of hot photoexcited carriers in these systems. We note that previous experiments on GaAs have observed the third 
order nonlinearity, but with an amplitude which is an order of magnitude smaller than predicted theoretically. 

For carbon nanotubes, one can anticipate at least three interesting applications of this phenomenon. First, as 
noted above, the method provides a means for current injection without electrical contacts. The absence of "low 
resistance" contacts on carbon nanotubes has often made it difficult to explore low energy transport phenomena in 
these systems |l4|-[l(| . A particularly interesting experiment would be to use the third order nonlinearity to produce a 
steady state separation of charge in a carbon nanotube rope or mat. In this state the "driving force" which produces 
a photocurrent via the third order nonlinearity would be balanced by the internal electric field produced by charge 
separation (in an open circuit condition). The relaxation of this charge distribution after the driving fields are turned 
off directly measures the conductivity along the pathways for charge motion in the system. Thus measurement of the 
transient relaxation after pulsed excitation would provide an interesting probe of the microscopic conductivity in this 
structurally heterogeneous system. Second, one can imagine applications which make use of the enhancement of the 
effect on semiconducting tubes (and its suppression in conducting tubes) to isolate semiconducting and conducting 
species in compositionally mixed samples. Finally, momentum transfer from the photoxcited carriers to intercalated 
ionic species can be used in principle to bias the diffusion of atomic or molecular species in the current carrying 
state. This effect requires in addition an asymmetry between the amplitudes for backscattering electrons and holes 
from the dopant species. This interesting application has been discussed in reference [pi using a simple model for the 
momentum transfer. 
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